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Abstract 

Angular X-ray cross-correlation analysis (XCCA) is an approach to study the structure of disor- 
dered systems using the results of x-ray scattering experiments. In this paper we summarize recent 
theoretical developments related to the Fourier analysis of the cross-correlation functions. Results 
of our simulations demonstrate the application of XCCA to two- and three-dimensional (2D and 
3D) disordered systems of particles. We show that the structure of a single particle can be recovered 
using x-ray data collected from a 2D disordered system of identical particles. We also demonstrate 
that valuable structural information about the local structure of 3D systems, inaccessible from a 
standard small-angle x-ray scattering experiment, can be resolved using XCCA. 
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I. INTRODUCTION 



Correlation methods are widely used in the physics of disordered materials such as 
amorphous and glassy systems [U |2]. Recently it was suggested [3H9] that angular cross- 
correlations of scattered intensities in coherent diffraction experiment could provide essential 
structural information about the system that is not accessible by other scattering techniques. 
The growing interest in the application of these methods is driven by two main factors. The 
first is its apphcation to the physics of disordered or partially ordered systems, and the sec- 
ond is an attempt to solve the structure of biological molecules in their native environment. 
Although ideas of applying cross-correlation techniques to study disordered systems go back 
to the work of Kam over 30 years ago (TUl E] , new perspectives are opened by the emerging 
high power free-electron lasers (FELs) which provide intense coherent femtosecond 

pulses of x-rays. 

In study of disordered systems application of cross-correlation methods is driven by the 
hope to unveil the hidden symmetries and dynamics in a disordered collection of elements 
composing the system. This leads to the problem of understanding systems with complicated 
correlation behavior between dynamical heterogeneity and medium-range order [T6HT8] in 
a large class of glass-forming liquids [T9142T] . One of the mysterious liquid systems whose 
structure is not fully understood is water. It exists in different forms known as low-density 
and high-density amorphous states [22] • The phase diagram of different forms is still under 
debate [23], El]. In such systems the relevant task could be, for example, the detection of 
an n-fold symmetry axis of an individual molecular species in a liquid composed of such 
molecules. An identification of bond angles or local structure in an amorphous system, from 
the study of angular correlations of the diffracted intensity is another attractive topic of 
research There are partly affirmative answers in the study of partially ordered quasi two- 
dimensional systems like liquid crystals [2S] in which hexatic bond order can be detected from 
the angular intensity correlation. Angular correlations with pronounced periodic character 
that can be directly related to the the local symmetry of the system were also observed in 
colloidal suspensions [3] |26| [27]. 

For biological systems it was proposed to use high power FEL pulses in so-called single 
particle coherent diffraction imaging experiments [28l |29] . In these experiments individual 
particles, for example macromolecules, are injected in the FEL beam and their diffraction 
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patterns are measured in the far-field. After the first experiments performed at LCLS in 
Stanford [3DJ, it was reahzed that even high power FELs provide single pulse resolution only 
up to tens of nanometers for gigantic viruses. Extension of the technique to sub- nanome- 
ter resolution for biological systems with the size below 100 nm is still under development. 
Clearly, it will be even more challenging to resolve molecules in their native liquid envi- 
ronment. One of the possible ways to enhance the scattered signal is the illumination of a 
large number of particles by a single PEL pulse. An approach based on cross-correlation 
analysis could provide a tool to determine three-dimensional (3D) diffraction pattern to high 
resolution with the available PEL sources [2IH33]- 

Here, we discuss potentials of application of x-ray cross-correlation analysis (XCCA) to 
the study of a certain class of disordered systems. Our model sample consists of identical 
non-interacting particles with orientational and positional disorder (see Pig. [T]). A coherent 
x-ray beam with a pulse duration shorter than the translational and rotational relaxation 
times of the system is scattered from the sample, and the diffraction pattern is measured 
in the far-field. The main concept of our approach based on application of XCCA to the 
study of disordered systems is depicted in Pig. |2j A large number of diffraction patterns is 
measured by successive pulses of the PEL beam with different realizations of the system. 
We assume that the time between the pulses is sufficient that two sequential realizations 
of the system are completely uncorrelated. XCCA is performed for each diffraction pattern 
in a sequence, and angular cross-correlation functions (CCPs) are averaged over the whole 
data set. In this way structural information about an individual particle in the system can 
be obtained. 

In this paper we summarize our recent theoretical results |lH6l [3l] based on the analysis 
of the high-order cross-correlation functions. It is organized as follows, in the next section 
the definition of the two-point and three-point CCPs, as well as their angular Pourier de- 
composition, is given. In the third section intensities scattered from our model system are 
defined. In section four a detailed description of XCCA for 2D dilute and dense systems 
is provided, this is followed in the next section by the description of a more general case 
of scattering from 3D systems. In section six major results of the recovery of the structure 
of an individual 2D cluster and analysis of the scattering from a 3D system of pentameric 
structures are discussed. This is followed by the conclusions and an overview section. 
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II. DEFINITION OF THE TWO- AND THREE-POINT CCFs AND THEIR AN- 
GULAR FOURIER DECOMPOSITION 



The two-point CCF defined for a single realization of a disordered system at two resolution 
rings qi and q2 is given by [U [TOl |35] [see Fig. ^c)] 

q2, A)= (7(gi, y.)J(g2, ^ + A)) , (1) 



where < A < 27r is the angular coordinate, I{q,(p) = I^q.'-p) ~ m.l^v))^^ is the intensity 
fluctuation function, {f{(p))^p = (l/27r) J^"^ f{ip)dip denotes the average over the angle ip. 

In a similar way, the three-point CCF for a single realization of a system is defined at 
three resolution rings gi, q2 and qs as [6] [see Fig. [3| 



C(gi,g2,g3,Ai,A2) = (/(gi,</^)/(g2,</' + Ai)/(g3,¥^ + A2)) , (2) 



where < Ai < 27r and < A2 < 27r are the angular coordinates. 

In practical applications one would need to consider CCFs {C{qi,q2, A))m and 
(C(gi, g2, ?3, Ai, A2))a/ averaged over a sufficiently large number M of diffraction patterns 
0, 

1 ^ 

(C(gi, g2, A))m = ^ 92, A)}"^, (3a) 

m=l 

(C(gi, g2, gs, Ai, A2))m = ^ 5^{C(gi, g2, ga, Ai, A2)}'", (3b) 

m=l 

where {C(gi, ga, A)}™ and {C(gi, ^2, ^3, Ai, Ag)}"" are the CCFs defined by Eqs. Q and ^ 
for the m-th realization of a disordered system. The importance of this statistical averaging 
will be discussed in the following sections of the paper. 

It is convenient to analyze the two-point CCF C{qi,q2, A) using a Fourier series decom- 
position in the (0, 27r) interval jl] 

00 

C{quq2,A)= J2 (4a) 

n=— 00 

Cai,.. = 7^ r C{qi,q2, A)e-^"^dA. (4b) 



27r 







Here, C^^ g^ is the n-th component in the Fourier series expansion of C{qi, q2, A), and Cg_^ g^ = 
at ?7. = by definition [see Eq. ([T])]. Substituting Eq. Q into Eq. (4b) and applying the 



Fourier convolution theorem we get 



/^n jn* jn (r\ 

91,92 - -'gi ■ -'92- W 

Here 1^ are the components of the Fourier expansion of the scattered intensity /(g, if) on 
the ring of radius q [see Fig. [3|^a)], 

oo 



/;^ = i-/ I{q,ip)e-^-^d^. (6b) 



n=— oo 
27r 







Since scattered intensities are always real quantities, it is easy to show that " = !„* and 



91:92 ~ 9l>92' 



In the specific case, when qi = q2 = q, Eqs. (4a) and ^ reduce to 



oo 



C(g,A) = 2j]C,"cos(nA), (7a) 

n=l 

c^=\i^\', c;>o. (7b) 



According to (7a) a strong single cosine dependence of C(g, A) can be observed for those 
values of q, at which one of the Fourier components significantly dominates over all 
others |3]. Such components can be related to the structure and symmetry of the system 

The Fourier series expansion of the three-point CCF C(gi, g2, 93, Ai, A2) can be written 

as 



00 00 



C{q„q„q„A„A,)= E ^9^93^'"'^'^'"'^'' (8^) 



Cl\Z, = (^] I I C^(gi,g2,g3,Ai,A2)e-"^^^e-™^^^ciAiciA2, (8b) 



ni=— 00 n2=—oo 
2 r2TT r2TT 



'0 JO 

where Cg^'^^^g^ are the Fourier components of the three-point CCF, and C^^'g^^g^ = for 
ni = 0, ^2 = and rii = —n2- Substituting Eq. ^ into Eq. (8b) one can get [6] 



9l>92,93 91 92 93 ' ^ / 

In general, Eq. ^ determines a relation between three different Fourier components of 
intensity of the order rii, 77-2 and ni + 77-2, defined on three resolution rings, gi, 52 and q^. 



In practical applications one would need to consider the Fourier spectra ((^5^,92 )m and 
^^qim,q3^^^ averaged over a large number M of diffraction patterns ^ E], 



M 



m=l 
M 



:iOa) 



(10b) 



m=l 



where {C^.^gJ"" and {C^ljg^'^^gJ"' are the Fourier components of the CCFs {C(gi, gs, A)}™ 
and {C(gi, g2, Ai, A2)}'" respectively, defined for the m-th realization of the system [see 



Eqs. (3a) and (3b) 



III. SCATTERING FROM A DISORDERED SYSTEM OF IDENTICAL PARTI- 
CLES 

We consider a scattering experiment in transmission geometry as shown in Fig. [l](a). A 
coherent x-ray beam scatters from a disordered sample, and a speckle pattern is measured on 
the detector in the far-field regime. In our simulations we consider a kinematical scattering 
approximation. As a general model system we assume a 3D sample consisting of identical 
particles with random positions and orientations. This model includes a variety of systems, 
for example, clusters or molecules in the gas phase, local structures formed in colloidal 
systems, viruses or complex biological systems in solution, etc. 

The amplitude Ak{q) scattered from the fc-th particle at the momentum transfer vector 
q can be defined as [1], 

^A:(q) = / Pfc(r)e*''-'-rfr, (11) 



where Pk{r) is an electron density of the k-th particle at the position [see Fig. [gb)] and 



the integration is performed over the volume of the particle. Using Eq. (11) the intensity 



/(q) coherently scattered from a disordered sample consisting of particles is given by 



N 



/(q)= e^^-'''=-'=^^l(q)^'^.(q) 



fcl,fc2 = l 

N 



= E / [ plM)Pk,{r,)e'''-'''>='-Mdr,dr,, (12) 

fcl,fe2 = l 

where the double summation is performed over all N particles, and the integration is per- 
formed over the volume of the /cj-th particle {i = 1, 2). Here, the following notation for the 



radius vectors is used, R-fcj.fci ~ ^k2,ki + ^21, where Rfc2,fei = R.fc2 — R-fci is the vector connect- 
ing two different particles ki and k2, and r2i = r2 — ri, where the vectors ri and r2 define 
the positions of scatterers (for example, colloidal spheres or atoms) inside the particles ki 
and k2 respectively [see Fig. [T]^b),(c)]. 

In the case of a partially coherent illumination and a dilute disordered system when the 
mean distance between the particles is larger than the coherence length of the incoming 
beam, the inter-particle correlations due to coherent interference of scattered amplitudes 



from the individual particles in Eq. (12) can be neglected. In these conditions, the total 



scattered intensity /(q) can be represented as a sum of intensities /^(q) = |v4fc(q)|^ corre- 
sponding to individual particles in the system 



TV 



/(q) = 5^4(q). (13) 



k=l 

IV. 2D DISORDERED SYSTEMS, SMALL ANGLE SCATTERING 

In this section we consider a particular case of a 2D system in a small angle scattering 
geometry when all vectors are defined in a 2D plane [see Fig. [l|^c)]. It was shown [4J, that in 
this case only even {n = 21,1 = 1, 2, 3, ... ) Fourier components of the intensity can have 
non-zero values. 



A. Dilute systems 



The intensity /^o(q) scattered from a single particle in some reference orientation ipo 
is related to the electron density of the particle p,/,o(i") through its scattered amplitude 



[Eq. (11)] as /^o(q) = I^Vo(^)l^- Similar to /(q) [see Eqs. (6a) and (6b)], the intensity 
J^(,(q) = /^g(g,(y9) can be represented as an angular Fourier series expansion. 



(14) 



where J^^^^ are the Fourier components of I^^^q^ip). 

For a dilute 2D system of identical particles, the intensity I^^,{q,ip) scattered from a 
particle in an arbitrary orientation ipk is related to the intensity I^g {q, (p) scattered from 
a particle in the reference orientation ipQ as I^^{q,ip) = /^^(g, — ipk)- Applying the shift 
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theorem for the Fourier transforms ^7] we obtain for the corresponding Fourier components 
of the intensities, = 1^^^ expl—inipk). Using these relations we can write for the Fourier 



components /" of the intensity /(g, (p) scattered from N particles [Eq. (13) 



9 

N 



jn _ jn p-inipk — jn \ (-1 r\ 



k=l 



where A„ = Ylk=i ^^v{~'^^'^k) is a random phasor sum 

The Fourier component /^^^ is related to the projected electron density py,o('^) ^f the 
particle as [H E] 

Oo = / y"p;„(ri)p^o(r2)Jn(g|r2i|)e-'"^-dridr2, (16) 

where is the angle of the vector Y21 in the detector plane, J„(p) is the Bessel function of 
the first kind of integer order n, and the integration is performed over the area of a particle. 
According to the structure of /^^^ its value strongly depends on the symmetry of a particle 
and determines selection rules ^ Ej for the values n of non-zero Fourier components ■ 
These selection rules can be used for the identification of the symmetry of particles in dilute 
systems. For example, for a cluster with 5-fold symmetry only n = 101, {I = 1, 2...) will give 
a non-zero contribution to the Fourier components of CCFs. 

Substituting Eq. (15) in Eq. ([s]), in the limit of dilute systems we have for the Fourier 
components C^^ ^^ of the CCF the following expression 



(-in _ rn* rn |A M fl 7^ 

The statistical behavior of A„ has been analyzed for different angular distributions of 
orientations of particles in the system [see Refs. [HE]]- It is clear, that in the case of a 
completely oriented system of particles (all tpk = 0) the square amplitude of the random 
phasor sum |A„| is equal to and C^^ ^^ = In the case of a uniform distri- 

bution of orientations of particles |A„|^ fluctuates around its mean value (|A,„|^) = N with 
the standard deviation cr|A„|2 = Averaging the Fourier components C^^ over a large 
number M of diffraction patterns decreases these fluctuations and leads to the following 
asymptotic result [U [S] 

( ,92 ) M = ^9i%0 C^o ■ ( I M ^-^9 1^0 C ,Vo • 0-^) 
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where (. . . denotes statistical averaging over M diffraction patterns. Importantly, the 
ensemble-averaged Fourier components {Cq^^q^ j^.^ converge to a scaled product of the two 
Fourier components of intensity Iq*^^^ and -^^j^^^ associated with a single particle. 

The amplitudes |/^^(,| and phases 0^ ,^^^ (for n 7^ 0) of the Fourier components /^^^^ = 



l-^g^Vol ^^P(^*^g.V'o) associated with a single particle can be determined using Eq. (18) [6]. This 
equation determines the phase difference between two Fourier components Iq^^^^ and Iq^^^^ 
of the same order n, defined at two different resolution rings qi and 

^rgKQ...>M]=€,^o-<.^o- (19) 

Similar to the Fourier components of the two-point CCF, the Fourier components of the 
three-point CCF Eq. (g can be expressed in the limit of a dilute system as |B] 

where A„j^„2 = exp{i[(ni-|-ri2)V'j~'^i'^j~'"'2'i/'fc]} is another random phasor sum. Our 

analysis shows [B] that in the case of a uniform distribution of orientations of particles the 
statistical average (A„^^„2)m converges to for a sufficiently large number M of diffraction 
patterns 

/rmi,n2 \ v t\t T{ni+n2)* jni Tn2 fn']) 



An important result of Eq. (21 ) is that the ensemble-averaged Fourier components {C^^^q^q^ 



91,92,93 /M 

converge to a scaled product of three Fourier components of intensity -^g"^^"^'**; -^gaVo' ^^'^ 
■^93^ ■00 dissociated with a single particle. Eq. (21) also provides the following phase relation, 



arg[(Q^,92:93>M] = €,.0 + €,^0 - (22) 

This equation determines the phase difference between three Fourier components 

Iql^^i and -^^3^^(j of different order n defined on three resolution rings. If rii = n2 = n 



and ns = 2n, equation (22) reduces to a particular form, giving the phase relation between 



Fourier components of only two different orders n and 2n. Phase relations (19) and (22) can 
be used to determine the phases of the complex Fourier components Iq^^^ using measured 
CCFs from a disordered system of particles [6]. 



B. Dense systems 



In the case of a dense system, when the average distance between particles is of the order 



of the size of a single cluster, the Fourier components of the intensity /(q) [Eq. (12)] can 
contain a substantial inter-particle contribution. In this case 1^ can be presented as a sum 
of two terms, 

/" = /part(?)+C-part(9), (23) 

where /partis') is attributed to a single particle structure discussed above, and /int-part(?) is 
defined by the inter-particle correlations [HE]. In a 2D system these two terms are [H [5] 

= I^^o^n, (24a) 
/rnt-part(?) = 2 5^ j J K,{r^)pk.Ar2)Jn{q\Rlik^)e "^^^'"^ dr.dr^, (24b) 

k2>ki 

where 0r2i ^ is the angle of the vector ^^^11^ the sample plane. 



Taking into account both terms of Eq. (23), the Fourier components of the two-point 



CCF for a single realization of a system ([s]) can be written as the following sum of four 
terms 

C^,,. = C ■ = Si + (^2 + ^3 ) + S2, (25) 

where 

5r=Ct(Si)/part(?2), (26a) 

^2 + ^3"= + iZ-^^AnKM. (26b) 

5r=/rnt-part(gi)/^t-part(?2). (26c) 

In Eq. (25) the term 5" is defined only by the structure of a particle, whereas the terms 



S'2 , 5*3 and 5*4 contain contributions from the inter-particle correlations due to the second 



term in Eq. (23). A similar expansion can be performed for the Fourier components of 
the three-point CCF. The influence of inter-particle correlations on CCF's measured on the 
same resolution ring were discussed in detail in Ref. [5]. 

V. 3D SYSTEMS, WIDE ANGLE SCATTERING 

In our previous discussion of x-ray scattering on 2D systems, we have seen that only even 
Fourier components of the CCFs have non-zero values. In the case of 3D systems non-zero 
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odd Fourier components can be also present when scattering to high angles is considered, 
due to Ewald sphere curvature effects. 

In general, the scattering vector q = (q"'",g^) can be decomposed into two components: 
q-*" that is perpendicular, and that is parallel to the direction of the incident beam [see 
Fig. |3](b)]. We also define the perpendicular R^^^j = 1^4,^1 + the z-components 

^kl,ki — ^k2,ki + Z21 of the radius vectors (see Figs. [T|a) and|3]). Using these notations we 



can write Eq. (12) in the following form 

TV 



/(q)= E e-^^^-^'=-'=i J y"p:^(r^,g^)p.,(4,g^)e^''"-^^."i^ir^^i4. (27) 

fcl,fc2 = l 

Here we introduced a modified complex valued electron density function, defined as 

Mri,qn = I PkArt,^)e-'''^'dz. (28) 



In the case of wide angle scattering, the effect of the Ewald sphere curvature [see Fig. [3| 
which manifests itself by the presence of the exponential factors e~*''^'^''2''=i and e"*''^^ in 



Eqs. (27) and (28), may become important. This effect can break the scattering symmetry 
of a diffraction pattern, characteristic for the scattering on a positive valued electron density 
(Friedel's law) and may reveal additional symmetries that can be hidden in the small angle 
scattering case. A wide angle scattering geometry may become important for scattering on 
atomic systems with local interatomic distances of the order of few Angstroms. 

For simplicity we will consider here a 3D system consisting of particles composed of 



identical scatterers. The modified electron density of a particle according to Eq. (28) can 
be defined in the following form 

Pk{^^, qn = f{q) E '^(^^ - (29) 

i=l 

where /(g) is a form-factor of a scatterer, and A^^ is a number of scatterers in the cluster. 
The coordinates (r^, Zi) define the position of the z-th scatterer inside the fc-th cluster. Using 



this definition and performing a Fourier transformation of Eq. (27) we obtain |3] 

N Ns _ 
ki ,k2 = l l,m=l 

where the summation over index / is performed over the positions of scatterers in the cluster 
ki, and the summation over index m is performed over the positions of scatterers in the 
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cluster ^2- We note here that due to the property of the Bessel functions [J«(0) = for 
n ^ 0] the terms with ki = ^2 and I = m are equal to zero. Taking into account that the 
terms with interchanged indices, i.e. ki,k2 and k2,ki, as well as l,m and m,l, differ from 
each other by a change of the sign of Z^J'i^_^ and by an additional factor (—1)", which arises 



due to the change of the phase 0R,xmi = ^r^'™ + tt, we have for even values of n in Eq. (30 ) 
i 

r(g^g^)=2(.n/(g)H Yl E ^^s J„(|q^| ■ |RSj)e" 

l<A:i<Ar l<l<Ns 
ki<k2<N l<m<Ns 

+ E E cos{q^Zf)U\q^\.\Ri-^\)e-''''^n (31) 



l<fe<Af l<l<Ns 
l<m<Ns 



and for odd values of n: 



r(g^g^) = -2(^rM/(?)l'[ E E ^in (g^Z™;,,) J„(|q^| ■ IRSD^" ' " 



k2,ki 



l<ki<N l<l<Ns 
ki<k2<N l<m<Ns 



+ E E sin(g^ZrVn(|q^|-|R^'|)e"^"'^^'"'], (32) 



l<A:<Ar l<l<Ns 
l<m<Ns 



where R^™' = R^fc^' and = Z^^. From the performed analysis we can see that, due 
to the curvature of the Ewald sphere (non-zero component), we obtain non-zero odd 
Fourier components of the CCF when scattering from a 3D system. These components 
become negligibly small for experimental conditions corresponding to a flat Ewald sphere, 
considered in the previous section. 



In the case of a dilute sample, the Fourier components of intensity defined in Eq. (29) 
reduce to 

N Ns 

nq\qn = {^T \f{<l)?Y. E e-^^^^^'jn(|q^| ■ |r^"'|)e-^"^^^-\ (33) 

k=l Lm=l 



VI. RESULTS AND DISCUSSION 

In this section we provide two examples of application of XCCA for investigation of the 
structural properties of disordered systems. 
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A. Recovery of the structure of an individual 2D particle 



In the pioneering work of Kam [TUl E] it was proposed to determine the structure of a 
single particle using scattered intensity from many identical particles in solution. However, 
this approach, based on a spherical harmonics expansion of the scattered amplitudes, was 
not fully explored until now. Recently, it has been revised theoretically [TJ [321 [35] and experi- 
mentally |8] . The possibility to recover the structure of individual particle was demonstrated 
in systems in two dimensions [TJ |Hl and three dimensions [32] using additional a priori 
knowledge on the symmetry of the particles. Unfortunately, these approaches, based on op- 
timization routines [32l [35] and iterative techniques [7l[H], do not guarantee the uniqueness 
of the recovered structure unless strong constraints are applied. Here we provide a brief 
overview of our recently developed approach [6J that enables a direct reconstruction of the 
single-particle intensity distribution using an algebraic formalism of two- and three-point 
CCFs without additional constraints. This approach is developed for a 2D system of parti- 
cles, but can be also used to study 3D systems provided that particles can be aligned with 
respect to a certain axis. 

Our goal here is to determine the scattering pattern of a single particle /^o(^) using a 



large number of diffraction patterns /(q) [Eq. (12)] corresponding to different realizations 
of the system. Once the single-particle intensity is obtained, conventional phase retrieval 
algorithms [39l [ID] can recover the projected electron density of a single particle. According 



to Eq. (14) the intensity scattered from a single particle can be uniquely determined by 
the set of complex coefficients {Z^^^} = {l-^^^ol' 'i^q.^/'o ~ ^^&i^q,ii>o)} ■ ^^^^ determine the 
Fourier components {Iq^^^} applying two- and three-point CCFs to the measured intensities 
I{q, if) scattered from particles. 



Equations (18), (19) and (22) constitute the core of our approach [6] and allow us to 
directly determine the complex Fourier components /^^g- According to our simulations, 
these expressions, originally derived in the dilute limit approximation [see the incoherent 



sum Eq. (13)], can be also used in the case of coherent scattering from a system of particles. 



The amplitudes l^g^^gl can be determined jS] using Eq. (18), and the phases (p^^^p^ using 
Eqs. (19) and (22), in particular when rii = n2. It might be tempting to use Eq. (7b) to 
determine the amplitudes by just taking the square root of (C^ )^^. However, in the 

case of coherent scattering from N particles the spectrum (C^)^ may contain a substantial 
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inter-particle contribution or noise [HE]. In contrast to this, even in the case of coherent 
scattering from a dilute system of particles, the inter-particle contribution to the spectrum 
i^qi q2) M determined at different resolution rings qi 7^ q2 is negligibly small. For the same 
reason we used the three-point CCFs [Eq. (|2])] defined at different resolution rings gi, q2 
and gs. Once all amplitudes and phases of the Fourier components /^^^^ are determined, the 



diffraction pattern corresponding to a single particle can be recovered using Eq. (14). 

We demonstrate our approach for the case of a coherent illumination of a disordered 
system of particles [Eq. (12)] with an incident fluence of 4 ■ 10^^ photons/ ^w? (that corre- 



sponds to 10^'^ photons focused on a sample area of 5 x 5 yum^), in the presence of Poisson 
noise in the scattered signal. As an example, we recover a diffraction pattern and projected 
electron density of an asymmetric cluster [Fig. |4]^a)] composed of PMMA spheres of 50 nm 
radius, with a size of a cluster d = 300 nm. We consider coherent scattering of x-rays with 
wavelength X = 1 A from a system of = 10 clusters in random positions and orientations, 
distributed within the sample area [see Fig. [l](c)]. Diffraction patterns are simulated for a 
2D detector of size D = 24 mm (with pixel size p = 80 /im), positioned in transmission 
geometry at L = 3 m distance from the sample [see Fig. [l|a)]. This experimental geometry 
corresponds to scattering to a maximum resolution of 0.25 nm~^. For given experimental 
conditions the speckle size corresponding to the illuminated area is below the pixel size of 
the detector. At the same time the speckle size corresponding to the size of a single particle 
is about 12 pixels, which provides sufficient sampling for the phase retrieval algorithm. 

The coherently scattered intensity simulated for a single realization of the system [41j is 
shown in Fig. 4]^b). The Fourier components of two-point and three-point CCFs, {C^^ 92 )m 
and (C*gi'^2,(?3 )j\f' were averaged over M = 10^ diffraction patterns. It is important to 
perform the averaging of the two- and three-point CCFs over a sufficient number M of 
diffraction patterns to achieve their convergence to the mean value. In practical applications 
the convergence of the CCFs can be directly controlled as a function of the number of the 
diffraction patterns M considered in the averaging [5J. The results presented in Fig. [5] 
demonstrate the evolution of the Fourier components {Cq^q^j^,^ as a function of M. As 
one can see, after averaging the CCFs over a sufficiently large number M of diffraction 
patterns, the amplitudes and phases of {Cq^^q.^ j^j asymptotically converge to their mean 
values. Comparing figures [5]^a),(b) and figures [5|^d),(e) one can see, that the CCFs converge 
slower in the presence of Poisson noise, and one needs to average over more patterns. This 
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effect is more prominent for the three-point CCFs. In Fig. [5[g)-(i) the phases of the Fourier 
components (Cg^'g^^g.^) are presented for three different combinations of rii and n2. It is 
well seen from these figures, that at the incident fiuence of 4 ■ 10^^ photons/ /im^ the phases 
and amplitudes of the Fourier components are reaching their mean values at M = 10^. 

The diffraction pattern corresponding to a single particle recovered by our approach is 
presented in Fig. ^c). As one can see, it reproduces the diffraction pattern of an individual 
cluster shown in Fig. |4]^a) very well. The structure of the cluster reconstructed from this 
recovered diffraction pattern was obtained by a standard phase retrieval approach [39] and 
is presented in Fig. |4]^d). The comparison of the single cluster structure obtained by our 
approach [Fig. |4]^d)] with the initial model shown in the inset of Fig. |4|^a) confirms the 
correctness of our reconstruction. The particle was reconstructed with a resolution of 25 nm. 
These results clearly demonstrate the ability of our approach to recover the single-particle 
structure from noisy data obtained in coherent x-ray scattering experiments. A detailed 
discussion of different aspects that can affect the quality of reconstruction, such as particle 
density, fluctuations of the number of particles N in the system, and noise can be found in 



B. Wide angle scattering from a 3D system of pentameric structures 

In the previous section we demonstrated the application of XCCA to the scattering 
data from a 2D disordered system of particles which led to a successful reconstruction 
of the structure of an individual particle. At the same time, for a 3D disordered system 
the same problem has been solved only if additional symmetry conditions were imposed 
[32|[33]. Hence, development of direct, self-consistent approaches for the recovery of the local 
structure of a 3D disordered system remains a challenge. As a step towards the solution of 
this problem we present the results of simulations of x-ray scattering experiments from 3D 
disordered systems composed of oxygen tetrahedral pentamers [Fig. [ij^b)]. The structure 
of an oxygen pentamer corresponds to the tetrahedral arrangement of water molecules in 
a so-called Walrafen pentamer jl3j, one of the favorable models of water that reproduces 
the results of x-ray and neutron scattering experiments P31 - H7] . We perform the XCCA of 
the diffraction patterns originating from such structures and show that the Fourier spectra 
of the CCFs contain information on the local structure of a disordered system that is not 
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accessible in the standard small-angle x-ray scattering (SAXS) experiment 

We consider here a scattering experiment on a 3D disordered system composed of 
oxygen pentamers [see Fig. [l|a)] with the following parameters: x-ray wavelength A = 1.0 A, 
detector size D = 24 mm (with a pixel size of p = 80 /im), sample-detector distance 
L = 25 mm. The samples were composed of a different number of oxygen pentamers. In 
our model we assume a uniform distribution of orientations of clusters. Three Euler angles 
were used to generate this distribution. The regular pentamer [see Fig. ^h)] was defined by 
the 0-0 distance of 2.82 A between the central oxygen atom and each atom in a vertex, and 
a tetrahedral angle 109.47° typical for water [H]. We demonstrate results of the calculations 
of the Fourier spectra (C^ and {Cq^^q^ averaged over a large number M of realizations 
of a system. All results presented in this section are normalized by the number of particles 
A^ considered in each particular case. 

First, we consider scattering from a dilute system of clusters, where the intensity /(q) 
scattered from A^ clusters is approximated by a sum of intensities scattered from each cluster 



[see Eqs. (13) and (33)]. The Fourier spectra (C^)^^^ defined on a single resolution ring q 
and averaged over M = 10^ diffraction patterns for A^ = 10 clusters in the system are 
presented in Figs. [6|^a)-(b). The spectrum C'^ for each realization of the system was directly 



determined using Eqs. (7b) and (33). Calculations were performed for a regular pentamer 
for the case of a flat (a) and curved (b) Ewald sphere. Calculations for the case of a flat 
Ewald sphere [IE] were performed by setting the z-component of the momentum transfer 



vector to = in Eq. ( 33 ) . In the case of a fiat Ewald sphere only Fourier components with 
even n- values have nonzero values. The amplitudes (C^ calculated for a single cluster 
A^ = 1 in the system (not presented here), and for a system containing A^ = 10 clusters 
[Figs. |6](a) and (b)], converge to the same functional dependence of q. Therefore, in dilute 
limit approximation the results do not depend on the number A^ of particles in the system, 
but are just scaled by A^. As one can see from Figs. |6](a)-(b), in the considered g-range 
only the Fourier components with n < 9 contribute to the spectrum, and in the case of the 
curved Ewald sphere the Fourier componets with odd n have magnitudes comparable with 
the ones with even n. Therefore, experimental conditions that correspond to curved Ewald 
sphere can provide additional structural information, as compared to the case with a fiat 
Ewald sphere. 

The effect of a distortion of the regular pentamer on the Fourier spectra of the two-point 
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CCF is demonstrated in Figs. [6[c),(d). These results were also obtained in the dilute limit 
using Eq. (33). The difference spectra (C^)^''^ 



' M 



dist 



Q I M 



reg 
M 



were calculated for 



the case when all atoms in each cluster were randomly displaced by 5% of the 0-0 distance. 
Here, {C^y^^ corresponds to the spectra calculated for the regular pentamer [Figs, [g] (a), 
(b)], and (C^)^^J* to the distorted one (not presented). These results show that CCFs 
preserve their functional dependence on momentum transfer vector q with small distortions 
of the cluster, and, importantly, the ratio between the values of different Fourier components 
is not significantly changed. 

It is interesting to consider similar results of simulations for the Fourier spectra {C^^ 



M 



7^a)-(d) the amplitudes {Cq^^q^) 



' M 



defined at different resolution rings qi ^ q2- In Figs, 
calculated as a function of q2 (qi = 1.71 A^^) are presented for a flat (a),(c) and curved 

'a),(b) were calculated in 



shown in Figs. 



(b),(d) Ewald sphere. The amplitudes {Cq^^q^) 
dilute limit from = 10 clusters, using Eq. (g in (g, and were averaged [Eq. ( |10a[ )] over 
M = 10^ realizations of the system. The spectra shown in Figs. [7]|^c),(d) were calculated for 
the case of coherent scattering from a cubic sample of 800 nm in size consisting of = 400 



clusters. In this case the coherently scattered intensities [Eq. (12)] were calculated for each 
realization of the system, and then the two-point CCFs and their Fourier components were 



determined using Eqs. (4a) and (4b) and averaged [Eq. (10a)] over M = 10^ diffraction 
patterns. A much larger number of diffraction patterns was considered in the latter case 
to average out the inter-particle contribution in the specra, which is absent in the former 
case. The difference between the spectra in Figs. [7|^c) and [Tj^a), and in Figs. [Tj^d) and 
^h), is presented in Figs. [7](e) and [7](f), respectively. It is easy to see from Fig. [7| that 
the inter-particle contribution to the Fourier components {C^^ 92)^/ ^^^^ coherent 

scattering from A^ = 400 clusters is negligible in the entire g-range, except of a small region 
corresponding to gi = g2 = 1-71 A~^. This is reflected in a sharp peak in the difference 
spectra presented in Figs. [7](e),(f). The width of this peak is related to the speckle size given 
by the sample size. As one can see from the insets in Figs, and[7](f), the magnitude of 
the inter-particle contribution at qi = q2 is much higher than the structural contribution. 
This clearly shows, that in the case of coherent scattering from a dilute system of particles, 
the inter-particle contribution to the Fourier components {Cq^^q^j^.j defined at qi 7^ q2 is 
negligibly small in contrast to the case qi = q2, where this contribution can be very large 
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Results presented here demonstrate the abihty of XCCA to obtain information about 
the oxygen clusters with tetrahedral arrangement of oxygen atoms typical for water. The 
spectra {Cg^ g^)j^ averaged for a large number of realizations of a system converge to a specific 
functional dependence of q, which one could hope to correlate with the local structure of 
the 3D system composed of many particles. Our results show, that small distortions of 
the oxygen pentamer do not lead to significant changes of the CCFs. This may help to 
identify an average local structure of water using experimentally determined CCFs. Our 
results demonstrate that the spectra {Cq^^q^)j^ defined at qi ^ q^ are more suitable for 
investigation of the structure of disordered system in the case of coherent scattering, than 
the spectra determined at q\ = q2. The Fourier components {Cl[ q.^) jyj can provide the 
structural information in the case when it is suppressed by a large inter-particle background 
in (^Cg)^. Our results show, that the Fourier components {Cg^^g^) can reveal structural 
information about local structure of a system with a larger number of particles, as compared 

VII. CONCLUSIONS 

In summary, we showed here how XCCA can be applied to study the local structure of 
disordered materials. Important steps in this analysis are Fourier decomposition of CCFs 
and averaging over many realizations of the system. As it was shown in our work, statistical 
fiuctuations of the Fourier components of CCFs, due to this averaging cancel out and infor- 
mation on the local structure of the system can be deduced. The most transparent results 
can be obtained in the analysis of the dilute 2D systems when contribution from the inter- 
particle correlations can be neglected. We demonstrated that in this specific case all Fourier 
coefficients of the angular decomposition of the coherently scattered intensity of a single unit 
can be completely recovered by the analysis of two- and three-point correlation functions. 
Importantly, not only amphtudes but also phases of this decomposition can be determined 
directly from the experimental data. Applying standard phase retrieval techniques to the 
recovered intensity gives the structure of a single unit in a disordered system. 

The XCCA of dense systems shows more comphcated features in the Fourier components 
of CCFs due to the presence of a substantial inter-particle contribution. It can be signifi- 
cantly reduced in some situations when CCFs are analyzed at different resolution rings. Our 
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analysis showed that inter-particle correlations may strongly influence the values of Fourier 
components of CCFs deflned at the same resolution ring. An extraction of a single-particle 
structure in dense systems by XCCA can be more complicated, however, it may provide 
information on the medium range order. 

The generalization of this approach to 3D systems is still a challenge. Here scattering to 
high angles can be especially interesting due to the presence of a signiflcant Ewald sphere 
curvature. Such scattering conditions allow to determine odd Fourier components of CCFs, 
which are inaccessible in a standard small-angle x-ray scattering experiment, and in this way 
explore hidden symmetries in the 3D system. This approach can be especially attractive for 
the study of atomic systems. 

We foresee that the general approaches described here for the analysis of the structure 
of disordered systems based on the angular cross-correlation techniques will find their wide 
applications in future. They could become especially attractive at a newly emerging free- 
electron lasers and diffraction limited ultimate storage rings |49j . 
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FIG. 1. (a) Geometry of the diffraction experiment. A coherent x-ray beam illuminates a 
disordered sample and produces a diffraction pattern on a detector. The direction of the incident 
beam is defined along the z axis of the coordinate system, (b) A disordered 3D sample composed 
of tetrahedral pentamers. All clusters have random position and orientation in the 3D space, (c) 
A disordered 2D sample composed of asymmetric clusters. All clusters have random position and 
orientation in the 2D plane. 
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FIG. 2. A concept of recovery of the structure of a single particle using x-ray scattering data from 
many particles. A large number M of realizations of a disordered system (a) composed of many 
identical particles is used to collect M diffraction patterns (b). X-ray cross-correlation analysis is 
applied to this x-ray data set to recover a diffraction pattern (c) corresponding to a single particle. 
The structure of a single particle (d) is determined applying phase retrieval algorithms to the 
recovered diffraction pattern (c). 
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FIG. 3. Scattering geometry in reciprocal space, (a) Scattered intensity I(q) defined in the 
detector plane in the polar coordinate system, q = {q,f)- (b) Ewald sphere construction. Here 
is the wavevector of the incident beam directed along the z axis, k°^* is the wavevector of the 
scattered wave with the scattering angle 2a. The scattering vector q = (q"*",*?^) is decomposed 
into two components that are perpendicular q-*" and parallel to the direction of the incident 
beam. Definition of the momentum transfer vectors in the derivation of the two-point C{qi,q2, A) 
(c), and three-point CCFs C(gi, 92, (Zs, Ai, A2) (d). 
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FIG. 4. (a) Scattered intensity (logarithmic scale) calculated for a single asymmetric cluster 
(shown in the inset), (b) Coherently scattered intensity from a disordered system consisting of 
= 10 clusters in random position and orientation. Scattered signal corresponds to the incident 
fluence of 4-10^^ photons//xm^ and contains Poisson noise, (c) Scattered intensity corresponding to 
an asymmetric cluster recovered from M = 10^ diffraction patterns of the form (b). (d) Structure 
of a single cluster reconstructed by an iterative phase retrieval algorithm using the diffraction 
pattern shown in (c). The intensities in (a), (c) are given in arbitrary units, and in (b) in photon 
counts. 
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FIG. 5. (a)-(f) Statistical convergence of the amplitudes 



c: 



11 ,12 I M 



and phases arg 



a 



11,12/ M 



at qi = 0.25 nm^^ q2 = 0.09 nm , n = 4,8, 14, and 16, calculated as a function of the number of 



diffraction patterns M used in the averaging, (g)-(i) The phases arg[(Cg/,q2^g3)jy^^] at qi = 0.23 nm~^, 
q2 = 0.24 nm~^ and q^ = 0.25 nm~^, calculated for three different combinations of ni and n2 as a 
function of M. Simulations were performed (a),(d),(g) without noise, (b),(e),(h) for the incident 
photon fluence 4 • 10^^ photons//xm^, and (c),(f),(i) for the photon fluence 4 • 10^^ photons///m^. 
In (b),(e),(h), and (c),(f),(i) Poisson noise was included. 
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FIG. 6. Fourier spectra (^C^'j j^^j calculated as a function of q for 1 < n < 9, for a flat (a),(c) and 
curved (b),(d) Ewald sphere. (a),(b) Amplitudes (C^ calculated in dilute limit approximation 
of scattering [Eq. (13)] from = 10 clusters, averaged over M = 10^ diffraction patterns. (c),(d) 



Difference spectra {Cg)'^^ calculated in the dilute limit approximation for a distorted oxygen 
pentamer for = 10, and M = (see text for details). All atoms in each pentamer cluster in 
the sample were randomly displaced by 5% of the 0-0 distance from the positions of a regular 
pentamer. 
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FIG. 7. Fourier spectra {Cq-^^q^) calculated as a function of q2 {qi = 1-71 A~^) for 1 < n < 9, for 
a flat (a),(c),(e) and curved (b),(d),(f) Ewald sphere. (a),(b) Amplitudes {Cq^^q,^) calculated in 
dilute limit approximation of scattering [see Eq. (13)] from = 10 clusters, averaged over M = 10^ 



diffraction patterns. (c),(d) Amplitudes 



C. 



91 .92 / M 



calculated for the case of coherent scattering 



[Eq. (12)] from a cubic sample of 800 nm in size consisting of A^ = 400 clusters, averaged over 



M = 10^ diffraction patterns. (e),(f) Difference spectra {Cq^^q^)'^^f show the difference (e) between 



the spectra in (c) and (a), and (f) between (d) and (b). The insets in (e) and (f) show full size of 
the peaks attributed to the inter-particle contribution. 
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